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ABSTRACT 



Context. Strong interactions between jets and stellar winds at binary system spatial scales could occur in high-mass microquasars. 
Aims. We study here, mainly from a dynamical but also a radiative point of view, the collision between a dense stellar wind and a 
mildly relativistic hydrodynamical jet of supersonic nature. 

Methods. We have performed numerical 2-dimensional simulations of jets, with cylindrical and planar (slab) symmetry, crossing the 
stellar wind material. From the results of the simulations, we derive estimates of the particle acceleration efficiency, using first order 
Fermi acceleration theory, and give some insight on the possible radiative outcomes. 

Results. We find that, during jet launching, the jet head generates a strong shock in the wind. During and after this process, strong 
recollimation shocks can occur due to the initial overpressure of the jet with its environment. The conditions in all these shocks are 
convenient to accelerate particles up to ~ TeV energies, which can lead to leptonic (synchrotron and inverse Compton) and hadronic 
(proton-proton) radiation. In principle, the cylindrical jet simulations show that the jet is stable, and can escape from the system 
even for relatively low power. However, when accounting for the wind ram pressure, the jet can be bent and disrupted for power 
:£ 10 36 ergs" 1 . 

Key words. X-rays: binaries - stars: individual: LS 5039 - Radiation mechanisms: non-thermal 



1. Introduction 

Strong shocks take place in very different scenarios like for ex- 
ample in extragalactic clusters of galaxies due to structure for- 
mation (e.g. Kang, Ryu & Jones 1996), active galactic nuclei 
(AGN) jets (e.g. Begelman et al. 119841 ), supernova remnants 
(e.g. van der Laan 1962), star forming regions (e.g. Bykov & 
Fleishman |1992| ), pulsar winds (e.g. Rees & Gunn |1974j l, galac- 
tic relativistic jets (e.g. Zealey, Dopita & Malin 1980), or young 
star jets (e.g. Araudo et al. 2007 and references therein). 

In the case of AGNs, the interaction of jets with their en- 
vironment allows the study of the jet properties, like its matter 
content and energy budget, either in the case of Fanaroff-Riley II 
(FRII; e.g. Scheck et al.Hg p Kra use[2005| and Fanaroff-Riley I 
(FRI; e.g. Perucho & Martf l2007i ) sources. The simulation of ex- 
tragalactic jets and their interaction with the environment have 
been addressed by a number of authors during the last decades 
(e.g., Marti et al. 119971 K omissarov & Falle 119981 Aloy et al. 
[19991 Rosen et al. [19991 Scheck et al. [20021 Leism ann e t al. 
120531 Perucho et al. [20(351 120061 and Mizuno et al. I20071) . In 
particular, the evolution of jets in homogeneous and inhomo- 
geneous ambient media, the influence of the internal energy on 
the structure of the jets, the mass load of external material, the 
evolution of FRI and FRII sources, and the growth of the Kelvin- 
Helmholtz instability in relativistic flows have been thoroughly 
studied. 

In the case of microquasars, X-ray binary systems with rel- 
ativistic jets (e.g. Mirabel & Rodriguez IT9991 Ribo 120051 ), there 
are not many studies on the consequences of the jet/medium in- 
teraction (e.g. Velazquez & Raga 2000). A common approach 
has been to compare with extragalactic jets to point out similari- 



ties between the two types of sources (e.g., Gomez 2001 , Hardee 
& Hughes 120031 ). It seems clear that these galactic outflows 
should also generate some kind of structure in their termination 
region, although the characteristics of this region from the ra- 
diative and dynamical point of view are unclear (e.g. Aharonian 
& Atovan [19981 He inz |2002l Heinz & Sunyaev EHOU Bosch- 
Ramon et al. 2005; Bordas et al., in preparation), and only a 
handful of sources show hints or evidences of jet/medium inter- 
action^. SS 433, Zealey et al. [19801 XTE J1550-564, Corbel 
et al. 120021 Cygnus X-3 , Hein dl et al. |2003l Cygnus X-l, Marti 
et al. 119961 Gallo et al. [20 051 H1 743-322, Corbel et al. l2005l 
Circin us X-l, Tudose et al. 120061 LS I +61 303, Paredes et al. 
120071 . 

A microquasar jet can interact strongly not only with the 
interstellar medium. In high-mass microquasars (HMMQ), the 
massive and hot primary star suffers severe mass loss in the 
form of a supersonic wind that can embed the jefl This wind 
can have strong impact on the jet dynamics at the binary system 
spatial scales. In fact, the jet may be destroyed by the wind, pre- 
venting jet detection. In such scenario, shocks would be likely 
to occur, leading to efficient particle acceleration and the pro- 
duction of multiwavelength radiation. All this shows the impor- 
tance of studying the interaction between the jet and the stellar 
wind using hydrodynamics. Previous works treated the jet/stellar 
wind interaction using a more phenomenological approach and 
focused on possible radiative outcomes (e.g. Romero et al. 2003 
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1 We focus on HMMQs due to the relevance of the stellar wind at 
the spatial scales studied here. The presence of dense material and the 
stellar photon field makes HMMQs even more interesting from the ra- 
diative point of view. In the case of low mass systems, there might be 
accretion disk winds affecting jet dynamics (e.g. Tsinganos et al. 2004), 
but the properties of the environment are not well known and the impact 
on the jet is less clear. 
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Romero et al. 2005 Aharonian et al. 2006a). Regarding the pro- 
duction of non-thermal particles in the present context, diffusive 
shock (first order Fermi) acceleration is a plausible mechanism 
(e.g. Drury 119831 1. Other mechanisms of particle acceleration 
may be operating as well, like second order Fermi acceleration 
(Fenni [19491 ), shear acceleration (e.g. Rieger & Duffy 2004]), or 
the so-called converter mechanism (e.g. Derishev et al. |20031 >. 

In this work, we present the results from 2-dimensional rel- 
ativistic hydrodynamical simulations, with cylindrical and slab 
symmetry, in a simplified scenario of the jet/stellar wind inter- 
action in an HMMQ. Our goal is mainly to make first order esti- 
mates of the dynamical impact of the wind on the jet. In addition 
to this, we provide with semi-analytical estimates on the radia- 
tive outcomes. 

In the next section, we briefly describe the main features of 
the scenario considered and present the physical parameters of 
the simulations; in Sect. [3] we characterize the setup of the simu- 
lations and explain their technical aspects; in Sect.|4] we present 
the results obtained, which are discussed in Sect. [5]together with 
the consequences derived from the simulations from the points 
of view of the radiation or jet stability. We summarize all this in 
Sect.|6] 

2. Physical scenario 

The physical scenario adopted here corresponds to an HMMQ 
with physical parameters similar to those of Cygnus X-l or 
LS 5039, two close high-mass X-ray binaries with jets and with 
moderate-to-strong stellar mass loss from a primary O star, in 
the form of a fast and supersonic wind. The distance between 
the jet base and the primary star is fixed to R or b = 3 x 10 12 cm. 
This distance is similar to those in LS 5039 (Casares et al. 2005 ) 
and Cygnus X-l (Gies & Bolton |T9861 >. Typical luminosities 
and photon energies for O stars are L„ « 10 39 erg s" 1 and 
eo ~ 3 kT « 10 eV, respectively; the stellar wind present speeds 
V w ~ 2 x 10 s cm s" 1 , mass loss rates M w * 10" 6 M yr" 1 , and 
temperatures T w ~ 10 4 K. 

For the jet, we adopt a hydrodynamical supersonic outflow 
ejected perpendicular to the orbital plane at mildly relativistic 
velocities, Vj.o ~ 10 10 cms _I , and temperature Tj a 10 9 " 10 K<k 
m p c 2 /kft, where m p is the proton mass, c is the speed of light 
and ks is the Boltzmann constant. The choice of the tempera- 
ture is somewhat arbitrary since our assumption is that the jet 
is supersonic. Nevertheless, this temperature plays a role and its 
relevance is discussed in Sect. In this work, we focus on jets 
formed by protons and electrons, as expected in case the jet is 
fed with particles from the accretion disk. Given the adopted 
temperatures and jet speeds, the most reasonable choice of adia- 
batic exponent is T = 5/3. 

If we considered a leptonic jet, the adopted temperatures 
would imply an important increase of the internal energy of the 
particles and thus would significantly reduce the Mach number, 
affecting the assumption of a supersonic flow used here. Scheck 
et al 120021 have shown that a supersonic leptonic jet with the 
same energy flux than a hadronic jet, develops in a very similar 
way. Therefore, large differences are not expected if the kinetic 
luminosity of the jet is kept constant regardless of the composi- 
tion. 

Regarding the geometry adopted for the simulation, a cylin- 
drical symmetry can suffice to study the jet head shock formation 
in the wind. Since the jet speed Vj » V w , the wind speed (which 
breaks the cylindrical symmetry) is neglected. As it will be seen 
in the slab simulations, this assumption is reasonable for power- 
ful jets, but not in the case of weaker ones. At the present stage, 



Table 1. Parameters of the wind 



Parameter 


Wind density (g cm" 3 ) 


2.8 x 10" 15 


Wind pressure (erg cm" 3 ) 


1.5 x 10" 3 


Wind velocity (cm s" 1 ) 


2 x 10 8 


Wind mass-loss rate (M Q yr" 1 ) 


10" 6 


Wind ram-pressure (erg cm" 3 ) 


1.12 x 10 2 



slab symmetry simulations are enough to show the importance of 
the stellar wind ram pressure on the jet, because the timescales 
on which the wind surrounds the jet (k. Z?j/V w ; where Rj is the 
jet radius) are longer than the simulation timescales. Also, the 
adopted configuration of the stellar wind, coming from a fixed 
side of the jet, is consistent with an orbital timescale much longer 
than the simulation timescales. The adopted simplifications are 
reasonable in our case, since we aim at pointing out remarkable 
features and not at obtaining accurate quantitative predictions. 

We note that the wind structure is a key point for our sim- 
ulations. It has been proposed that massive stellar winds may 
be porous (e.g. Owocki & Cohen 2006). In addition, mass-loss 
rates of OB stars may be smaller up to a factor of several due 
to wind clumping, as shown, e.g., in Puis et al. (2006), where 
estimates for the real mass-loss rates for several O stars after 
wind dumpiness correction are given. For very porous winds, 
the collision of the jet with very small and dense clumps (i.e. the 
dynamically relevant part of the wind) would be rare, although 
these interactions could have interesting radiative outcomes (e.g. 
Aharonian & Atoyan 1996; Romero et al. 2007a). Also, smaller 
mass loss rates would weaken to some extent the jet kinetic lu- 
minosity constraints derived from this work. 

In our simulations, the jet starts at a distance zo = 6xlO I0 cm 
from the compact object, where it should be still little affected 
dynamically by the environment. The initial half-opening angle 
of the jet is taken to be 0. 1 radians, which is given by the rela- 
tion between the advance and expansion velocitiefl Then, the 
jet propagates through the system sweeping up stellar wind ma- 
terial. Two shocks are eventually generated, one in the stellar 
wind (i.e. forward or bow shock) and another one in the jet itself 
(i.e. backwards or reverse shock). The shocked material is called 
here shell for the wind (i.e. bow shock downstream) and cocoon 
for the jet (i.e. reverse shock downstream). The jet continues 
its advance slowing down due to kinetic energy exchange with 
the swept up wind material, and eventually could be bent, and 
even destroyed, by the lateral stellar wind ram pressure. Also, 
the pressure of the wind or the cocoon as compared to the lateral 
pressure of the jet can give rise to strong recollimation shocks. 
Recollimation shocks may also destroy the jet via deceleration 
and loss of collimation of the fluid downstream of these shocks 
(Perucho & Martf 120071 . 

3. Jet/wind interaction simulations 

We have performed five numerical simulations, three with cylin- 
drical symmetry and two with slab symmetry, in order to study 
the evolution of HMMQ jets with different injection power. 

2 A jet in free expansion due to a large initial overpressure expands 
approximately at its sound speed (1991 1, which is, in our simulations, 
roughly an order of magnitude smaller than the advance speed (for the 
temperatures given in Sect.[3}, which results in the opening angle given 
above. 
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Table 2. Parameters of the jets 



Parameter 


I 


II 


III 


Jet power (erg s" 1 ) 


3.0 x 10 34 


10 36 


3.0 x 10 37 


Jet init. pressure (erg cm" 3 ) 


9.1 


68 


6.2 x 10 4 


Jet init. density (g cm" 3 ) 


2.2 x irr ls 


5.9 x 10" 16 


1.8 x 10" 14 


Jet init. temperature (K) 


7.4 x 10 8 


2.1 x 10 9 


6.2 x 10 10 


Jet init. speed (cm s" 1 ) 


1.3 x 10 10 


2.2 x 10 10 


2.2 xlO 10 . 


Jet init. Mach number 


47.02 


51.05 


9.35 



The simulations were performed using a 2-dimensional 
finite-difference code based on a high-resolution shock- 
capturing scheme which solves the equations of relativistic hy- 
drodynamics written in conservation form. This code is an up- 
grade of the code described in Marti et al. (1997) (e.g. see 
Perucho et al. 2005). Simulations were performed in two dual- 
core processors in the Max-Planck-Institut fur Radioastronomie. 

The numerical grid of the cylindrical geometry simulations 
is formed by 320 cells in the radial direction and 2400 cells in the 
axial direction in an uniform region, with physical dimensions of 
20x 300 Rj. An expanded grid with 160 cells in the transversal 
direction brings the boundary from 20 Rj to 60 Rj, whereas an 
extended grid in the axial direction composed by 480 extra cells 
spans the grid axially from 300 Rj to 450 Rj. The enlargement of 
the grid is done to take the boundary conditions far enough from 
the region of interest and avoid numerical reflection of waves in 
the boundaries affecting our results. The numerical resolution in 
the uniform grid is thus of 16 cells/Rj in the radial direction and 
8 cells//?] in the axial direction. Outflow boundary conditions 
are used on the outer boundaries of the grid, inflow at injection, 
and reflection at the jet axis in the cylindrical case. In the sim- 
ulations, all the physical variables are scaled to the units of the 
code, which are the jet radius, Rj, the rest-mass density of the 
ambient medium, and the speed of light. 

The jet is injected in the grid at a distance of 6 x 10 10 cm 
from the compact object, and its initial radius is taken to be 
^j,o = 6 x 10 9 cm. The time unit of the code is thus equivalent to 
0.2 seconds, as derived from the radius of the jet at injection and 
the speed of light (R^q/c). The grid covers the distance between 
6 x 10 10 cm and 2 x 10 12 cm, i.e., a 60% of R mb . The ambient 
medium, i.e. the stellar wind, is composed by a gas with thermo- 
dynamical properties derived from Sect. [2] (see Table [TJ. Both 
the jet and the ambient medium are considered to be formed by 
a non relativistic gas with adiabatic exponent F = 5/3. 

In the simulations the jets are injected with different prop- 
erties. The physical parameters that characterize each simula- 
tion are listed in Table [2] the main difference between the three 
cases being the kinetic luminosity (Lj): weak jet (case I) of 
3 x 10 34 erg s -1 ; mild jet (case II) of 10 36 erg s _1 ; powerful jet 
(case III) of 3 x 10 37 erg s" 1 . The selection of the jet power is 
addressed to show some illustrative cases: the lowest value, that 
of case I, is similar to the minimum power required to power ra- 
dio emission under reasonable assumptions of the radio emitter 
(e.g. for Cygnus X-l, see Heinz 2006); that of case II is a bit 
smaller than the power estimates found for Cygnus X-l (Gallo 
et al. 120051 ), and similar to the Lj inferred for LS 5039 (Paredes 
et al. 120061 ), and the maximum value, that of case III, is between 
the upper limit for the Cygnus X- l Lj an d the Lj of SS 433 (e.g. 
Gallo et al. 120051 Marshall et al. 2002). These jets are evolved 



until they have reached distances similar to the binary system 
size. 

The velocities of the jets were selected as those of mildly 
relativistic flows, with Lorentz factors y s = 1.1 for case I, and 
yj = 1.5 for cases II and III. The jet Mach numbers are similar 
in cases I (47.02) and II (51.05), whereas it is much smaller in 
case III (9.35). In the latter, the increase in the jet power is due 
to an increase in the internal energy of the jet, keeping the same 
jet speed as in case II. The jet temperatures also show this fact, 
as the jet in case III has a temperature larger than that of case II 
by a factor 30. This allows us to see the influence of changes of 
velocity and internal energy on the jet evolution. 

In the case of the slab geometry simulations, the grid size is 
shortened to 200 Rj along the jet axis, and it is doubled in the 
direction transversal to the jet axis (it extends from -60 Rj to 
60 Rj). The numerical resolution is the same as that used in the 
cylindrical simulations. The jet is injected in an ambient medium 
that mimics the presence of a spherical wind centered in a star 
at a distance R ot \, from the compact object on the orbital plane, 
with the properties mentioned in Sect. [2] In this case, the bound- 
ary conditions are injection on the side from which the wind is 
assumed to come and on the base of the jet, and outflow in the 
outer boundaries on the opposite side of the grid and at the end 
of the grid in the axial direction. This kind of simulations were 
done for cases I and ifl since the interaction with the wind was 
expected to be possibly relevant for the jet dynamics. 



4. Results 

4. 1. Case I: a weak jet 

Cylindrical jet 

The jet in case I needs ff ~ 900 s to cover the distance be- 
tween the injection and the outer boundary of the numerical grid. 
The bow shock generated by the injection of the jet in the am- 
bient medium moves at a mean speed of Vb s ~ 0.06 c. The jet 
generates a high pressure cocoon that keeps it collimated (see 
the discussion in the next section). No other shocks or signifi- 
cant jet perturbations are generated, and the jet propagates well 
collimated through the wind. We note that the speed of the un- 
shocked jet does not change significantly along the axis. We do 
not show any image of this simulation due to the lack of remark- 
able features. A similar result has been obtained in case II, we 
thus refer the reader to Fig. |2](see Sect. |4.2t . 

The approximation to a static wind is reviewed critically in 
view of the results of the slab jet simulation presented next, in 
which the wind has the appropriate velocity. 

Slab jet 

The same simulation has been performed in planar coordi- 
nates including the wind motion. The wind comes from a point 
source at the distance and position of the primary star. Figure [TJ 
shows panels of axial velocity, rest mass density and pressure 
for the moment when the jet head has left the grid. The impact 
of the wind in the evolution of weak jets is clear from the plots. 
The bow shock generated by the jet moves at 0.05 c and is de- 
formed by the wind and the asymmetry in pressure on both sides 
of the jet. The wind also pushes the jet sideways, and two diag- 
onal shocks at z ^ 3 x 10 cm are formed due to the interaction 
with the environment, which is in overpressure with respect to 
the jet flow. After the first shock, the jet widens and recollimates, 
starting a new process of expansion after z ~ 2.6 x 10 cm. The 



3 The slab jets have the same lateral pressure and Lj/cm -on the sim- 
ulation plane- as in the cylindrical case. 
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jet is finally disrupted due to the action of the wind, after the 
jet has lost some of its initial inertia due to the previous inter- 
nal shocks. A strong shock is generated in the disruption point 
at z ~ 5 x 10 n cm. Farther from this point the jet is signifi- 
cantly deviated and a portion of fast material of the remaining 
outflow is observed to leave the regular numerical grid sideways 
at z ~ 6 - 7 x 10 11 cm. This simulation shows an example of a 
frustrated jet turned into a decollimated, slow outflow due to the 
effect of the wind. 



the stellar wind, recollimation shocks may change their posi- 
tion and properties. However, as is shown in the next section, 
these shocks may still be present. Such shocks, if strong enough, 
can trigger pinching of the jet, as observed in Fig. [4] which, in 
turn, may end up in mass load of ambient gas and jet disruption 
(Perucho & Marti 2007 ). Hence the importance of these recolli- 
mation shocks, not only from the point of view of particle accel- 
eration and subsequent radiation, but also from the dynamics of 
the jet. 



4.2. Case II: a mild jet 



5. Discussion 



Cylindrical jet 

The simulation of case II in cylindrical coordinates is very 
similar to that of case I. Fig.|2]shows the rest-mass density and 
pressure of the jet when its head has reached the end of the grid. 
The main difference between this case and case I lies on the fact 
that here some internal shocks are produced by the slight over- 
pressure of the jet with respect to the cocoon at z > 1.4xl0 12 cm. 
Before, the jet is efficiently collimated by the high pressure co- 
coon. The bow-shock propagates at Vb. s ~ 0.3 c, so it is signifi- 
cantly faster than in the previous case. The bow-shock introduces 
a jump in pressure of 10 8 with respect to the cold wind. The in- 
ternal, conical shocks observed in the jet itself, show jumps in 
pressure of 10 2 - 10 4 , although there are no strong changes of 
axial velocity with respect to the unshocked jet material. This 
implies that not much kinetic energy could be transferred to non- 
thermal particles via these shocks. 

Slab jet 

Also in case II, a slab simulation was performed account- 
ing for the velocity of the stellar wind from the primary star. 
Fig.[3]shows several maps of physical magnitudes from the sim- 
ulation. The bow-shock is also faster than that in the slab case 
I (Vt>s ~ 0.3 c). Several diagonal shocks are formed due to the 
interaction with the wind along the jet trajectory. Again, the jet 
is in underpressure with respect to the ambient, what generates 
an initial collimation. The pinching triggered by the first collima- 
tion shock and subsequent expansion of the jet, together with the 
thrust exerted by the wind on the jet, generate nonlinear struc- 
tures on the latter. In this case, the jet is also finally disrupted at 
z ~ 10 12 cm (see Fig. [3]). The pattern of shocks at the disruption 
region is similar to that observed in Fig. Q] These results show 
that the jet can be disrupted by the stellar wind ram pressure even 
for intermediate kinetic luminosities. 

4.3. Case III: a powerful jet 

The simulation of the powerful jet (Lj = 3 x 10 37 ergs _1 ) 
shows important differences with respect to the previous cylin- 
drical simulations. Fig. [4] displays the maps of rest-mass den- 
sity, pressure and axial velocity when the bow-shock is about to 
reach the limit of the regular grid. The bow-shock propagates at 
Vbs ~ 0.3 c. This Vbs is the same as in case II because part of the 
larger initial thrust of case III jet is invested in lateral expansion 
and heating. The latter happens due to, in contrast to the previous 
cases, a strong recollimation shock formed at z ~ 1.3 x 10 12 cm. 
This shock is formed due to the initial overpressure of the jet 
with respect to the cocoon. 

The relation between the pressure in the cocoon and that in 
the jet can make possible the formation of such shocks. The jump 
in pressure generated in the jet axis is of the order of 10 4 . In 
a second stage of the jet evolution, i.e., once the bow-shock is 
away from the studied region and the cocoon is substituted by 



5.1. On the formation of shocks in the studied region 

We see that two different kinds of shock, important for the ac- 
celeration of particles and subsequent radiation, may form: re- 
verse/bow shocks, which cross the region of constant density 
(within ~ R OI b) at a fraction of c on timescales < 10 4 s; and rec- 
ollimation shocks, generated by the interaction of the jet with its 
environment, the cocoon or the stellar wind itself. Unlike the re- 
verse/bow shock, which appears naturally as shown by the sim- 
ulations, the formation of recollimation shocks depends on the 
complex jet/medium pressure balance and deserves further anal- 
ysis. The possibility that they arise is discussed in this section in 
terms of the properties of the jets and their environments. 

Recollimation shocks form when an initially overpressured 
jet expands and reaches a lower pressure than that of its envi- 
ronment. In order to know whether a jet generates such a shock, 
and at which distance, the change of jet pressure with distance is 
needed to be compared with that of the ambient medium. In an 
adiabatic expansion, the condition pjjjVjAj = C\ holds, where 
C\ is a constant and Aj is the jet cross section. In our simula- 
tions, the jet velocity and Lorentz factor are basically constant, 
implying pj oc l/Rr, which, for a conical jet, transforms into 

Pj oc 1/z 2 , i.e., the density of the jet decreases as the square of 
the distance to the source. It is also known that the jet pressure 
(Pj) and density satisfy the relation Pj/p^ = C2, where T is the 

adiabatic exponent and C2 a constant. Thus, Pj oc 1 /z 2r , which 
implies (F = 5/3): 



P oc z 



-10/3 



(1) 



A minimum distance for the position of a recollimation shock 
can be given as the place where the pressure in the jet is the 
same as the pressure of the ambient (z eq )-0We can express this 
distance in terms of the pressure of the jet at a reference position 
Zo (Pjfi), which will be taken as the injection point in the grid, 
and the pressure of the ambient medium (P ex t)- From Eq. ([TJ: 



Zs > 



J.<> 



3/10 



(2) 



From this, certain limits on the pressure (and temperature) of 
jets, suitable to develop recollimation shocks, can be given. 



4 Actually, the shock occurs at a distance given by z sq , the Mach 
number of the jet at this point (M), and the jet radius at z cq (R cq ): 
z ~ Zcq + ^eq x M (see 11991) , as the information takes a certain time 
to reach the jet axis. At this distance, the pressure in the jet is strictly 
smaller than that of the ambient. In fact, the distance at which the cur- 
vature of the jet starts to change, which can be identified with the place 
where P s ~ /> cxt , is about half of that in which the shock occurs (see 
Fig- EJ- The second part of the equation depends on the parameters of 
the jet in the equilibrium point. For simplicity, we adopt the z cq as a 
lower-limit approximation to the position of the shock. 
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Fig. 1. Maps of axial velocity (cm s" 1 ), rest mass density (g cm" 3 ) and pressure (dyn) at the end of the simulation (ff = 680 s) 
of a weak jet in slab geometry, case I. The horizontal and vertical coordinates indicate distances (cm) to the injection point in the 
numerical grid and the jet axis, respectively. The stellar wind is centered on the position of the primary star, at x — -6 x 10 10 cm 
and y = 3 x 10 12 cm in the figure. 



Since we are interested in the region of maximal wind/jet in- may result in a recollimation shock in < R OI \, can be given. If zo 

teraction, we impose that the right hand side of the inequality in is taken as the injection point, i.e., zo = 6 x 10 10 cm, we obtain: 

Eq. d21> is < 10 12 cm, being a safe estimate for z s < R or \,. Thus, p. 

fixing zo, an estimate of the initial overpressure in the jet that 1 < — — < 1.2 X 10 4 . (3) 

"ext 
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Fig. 2. Maps of rest mass density (g cm 3 ) and pressure (dyn) at the end of the simulation (ff = 164 s) of the mild jet (case II) in 
cylindric geometry. The horizontal and vertical coordinates, and the star location, are like in Fig.Q] 



Eq. (O can be written in terms of the jet temperature (7j o) and 
the jet density (pj,o) at the injection point: 



T i,o< 



1.2 x 10 4 f ext m p 

(T - l)pj,0&B 



K. 



Pi,o 



(4) 



Both cylindrical and slab simulations show that, initially, the 
jets are embedded in their cocoons. However, as shown by the 
slab jet simulations, quickly the stellar wind starts to dynam- 
ically affect the bow-shock and the jet itself, and once the jet 
head has left the binary system, shocked wind will substitute co- 
coon material. We note that the properties of the wind change 
via shocking with the jet. In this case, the relative pressure of 
the jet with respect to the ambient would also change, at least 
on the side where the wind impacts. This asymmetry in the am- 
bient medium can give rise to asymmetric shocks, generated on 
one side of the jet and propagating diagonally through it, as ob- 
served in the slab simulations (Figs.Q]and[3]). In short, these two 
different scenarios, i.e. interaction with the cocoon and the wind, 
have to be taken into account. 

When the pressure of the ambient medium corresponds to 
that of the cocoon, the following equation (Begelman & Cioffi 
1989) estimates the pressure of the latter (P c ) as a function of 
the jet kinetic power, advance speed of the bow shock, and the 
surface of interaction between the bow shock and the ambient 
medium (A c ): 



(5) 



V bs A c 



When the external pressure is exerted by the wind on the jet, it 
can be computed by assuming that the density of the shocked 
wind is four times the original density of the wind, basing on the 
Rankine-Hugoniot jump conditions, and that most of the kinetic 
energy of the wind is thermalized and, thus, its internal energy 
is similar to the original kinetic one. 
Jet/Cocoon interaction 

For a typical jet with kinetic power 10 36 ergs~', bow-shock 
advance speed 0.1 c and bow-shock radius around 3 x 10 n cm, 
the cocoon pressure is P c ~ 10 3 erg cirT 3 . Then, if the density in 
the jet at injection (zo — 6 x 10 10 cm) is pj o ~ 10~ 15 gcirr 3 , the 



temperature must be: 



7j, < 3 x 10 



14 



10 36 erg s 
3 x 10 n cm\ 2 



3 x 10 y cms 

10 _15 gcnT 3 

Pj,o 



K, 



(6) 



which is fulfilled accounting for the fact that the jet must be su- 
personic (i.e. Tj o < (F - l)m p c 2 /&B ~ 10 12 K). 

One could also impose that the pressure at zo must be larger 
than that of the ambient to have an initial expansion of the jet 
(see Eq. [3j. A lower limit in the temperature of the jet may be 
set: Tjfl > 7j ; o ; max/(l-2 x 10 4 )K. In the cylindrical simulations 
of cases I and II the temperatures of the jets are below this limit 
(see Table [2}, so they do not generate recollimation shocks (see 
Fig. Nevertheless, we note that the location of the injection 
point in the numerical grid, zo, is somewhat arbitrary. As the 
pressure in the jets increases towards the compact object like 
z 51 ^ 3 , it can be inferred that, with the physical conditions de- 
rived from the simulations, the pressure in the jet would, at some 
point, become larger than in the cocoon. Hence, the appropriate 
conditions for the generation of a recollimation shock in those 
simulations would then exist. 

Jet/wind interaction 

In the case of the jet interacting with the stellar wind, 
the pressure is derived as explained above to obtain P w ~ 
10 2 erg cirT 3 . Using the jet density given above for the jet, 



7j <0 < 3 x 10 1 



2.8 x 10- 15 gcirr 3 / \2xl0 8 cms 
lO^gcirT 3 

Pj.o 



4 



K, (7) 



which is again fulfilled accounting for the fact that the jet must 
be supersonic. 

We conclude that the presence of standing shocks within the 
binary system is very probable within our prediction limits. 
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Fig. 3. Maps of axial velocity (cm s _1 ), rest mass density (g cirT 3 ) and pressure (dyn) at the end of the simulation (tf = 140 s) of a 
mild jet in slab geometry, case II. The horizontal and vertical coordinates, and the star location, are like in Fig.Q] 



5.2. Radiative processes the time they remain within the system, since the shock energy 

transfer reduces further out with the strong decrease of the wind 
As discussed above, strong shocks take place in the jet head, i.e. density. In addition, quasi-permanent recollimation shocks be- 
the bow/reverse shocks, and in the jet/cocoon contact surface, tween the jet and the stellar wind can also occur, as noted in 
i.e. the recollimation shocks. In practice, the capability of the Sect. O Since all these shocks are strong, and we know ap- 
jet head structures to produce non-thermal emission is linked to proximately the speed of the shocks from the simulations, we 
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Fig. 4. Maps of axial velocity (cm/s), rest mass density (g cm 3 ) and pressure (dyn) at the end of the simulation (ff = 212 s), for a 
powerful jet (case III). The horizontal and vertical coordinates, and the star location, are like in Fig.Q] 



can estimate the particle acceleration efficiency. In the follow- 
ing, we compute the maximum energy and the emission spectral 
energy distribution of particles accelerated under shock condi- 
tions similar to those shown by the simulations. 

5.2.1. Non-thermal particle production 

Particle acceleration 

As a first order approach, and provided that the shock 
Lorentz factors are ~ 1, we compute the particle acceleration 
rate in the context of the non-relativistic (first order Fermi) dif- 
fusive shock acceleration mechanism, in the test particle ap- 
proximation, 1 -dimensional case and with dynamically negligi- 
ble magnetic field. This mechanism typically predicts, for strong 
shocks, a power-law particle distribution of index -2 (e.g. Drury 
1983 ). The acceleration rate in cgs units can be computed using: 



E = rjqBc, 



(8) 



where E is the particle energy gain per time unit, q the particle 
charge, B the magnetic field in the accelerator and rj the acceler- 
ation efficiency coefficient, which can be calculated through the 
formula (e.g. Protheroe 1999): 



>1 



(9) 



for a parallel magnetic field, where is proportional to the dif- 
fusion coefficient (K^ = 1 corresponds to the Bohm regime) and 
Vs is the speed of the shock in the upstream reference frame. The 
value of the second parameter can be obtained from the simula- 
tions, allowing to estimate rj. For the reverse and forward shocks, 
we have taken typical velocities when the jet head is still within 



the binary system, obtaining Vs ~ 10 ll) cm s . In case of recol- 
limation shocks, as shown in Fig.|U velocity jumps could reach 
similar values. In the Bohm diffusion case, r\ ~ 2 x 10~ 2 . We 
recall that the recollimation shock shown in Fig. 0] is related to 
the cocoon, although such structure could also be produced by 
interaction of the jet lateral surface with the stellar wind. 

Concerning the B strength, the hydrodynamical assumption 
imposes an energy density for the adopted magnetic field <K the 
(internal plus kinetic) matter energy density. Thus, the value of 
B is to be well below equipartition with jet matter to be consis- 
tent with our simulations. For the (internal plus kinetic) matter 
energy density in the shocked material given by the simulations, 
which ranges ~ 10 2 - 10 6 erg cm 3 , B should be <k 10 2 - 10 4 G 
depending on the jet power and shock properties. 

We can compute the maximum energy of the particles equat- 
ing Eq. (O to the particle energy losses. 

For electrons, when the synchrotron energy losses are domi- 
nant, the maximum energy can be obtained from: 



(10) 



For inverse Compton (IC) losses in the Thomson regime (E <sc 
5.1 x 10 s eV/eo; eo = 10 eV), and accounting for the stellar 
radiation density in the jet, U* ~ 300 erg ctrT 3 , the maximum 
particle energy is: 



V^TeV, 



(ID 



which is to be substituted by the following formula in the Klein- 
Nishina (KN) regime (E > 5.1 x 10 5 eV/eo) (Khangulyan, 
Aharonian & Bosch-Ramon 2008): 

£max « 10 9 (B77)" TeV. (12) 
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Zsmax is determined by the smallest among the shown E m:ix . 
For K& = 1,77 = 2xl0~ 2 , synchrotron energy losses are dominant 
for B > 0.2 G. £ max can reach 1 TeV for B = 50 G, and 10 TeV 
for B = 0.4 G. 

In the case of inefficient energy losses, as it commonly hap- 
pens for protons, the size of the accelerator, approximated here 
as the diameter of the jet at the shock location (Djs), gives the 
very upper-limit to the maximum particle energy when equated 
to the particle mean free path (the gyroradius in our case): 

£>j,s = r g = — (13) 
qB 

Higher energy particles escape the accelerator since they cannot 
be confined by the magnetic field. For Dj s ~ 10 cm, E m . dx may 
reach ~ 10 3 TeV. 

Characterizing the particle energy distribution 
Once electrons are injected in the emitter, they suffer losses 
via synchrotron and IC processes. Relativistic Bremsstrahlung 
and ionization losses of electrons are negligible due to the long 
timescales and the high ionization degree of the shocked mate- 
rial, respectively. Protons lose energy mainly via proton-proton 
collisions, although the efficiency of this process is in general 
small. All the particles are convected with timescales that are 
shorter than the radiative ones for the electrons with the higher 
energies, and protons of any energy. The convection timescale is 
approximated by r conv ~ ld/Vd ~ 400 s, being /<j ~ 10 12 cm and 
Vd ~ Vs/4 (from the strong shockjump conditions) the length of 
the downstream region and the convection velocity, respectively 
(derived from the simulations). 

In order to respect the restriction of a cold jet and test particle 
shocks, we assume that the non-thermal particle energy density 
is smaller than the thermal one, fixing the ratio between them 
to 0.1 very close to the shock surface. Accounting for this, and 
under the conditions of the shocks considered here, i.e., a power 
similar to that of the jet going to thermal particles downstream, 
the injected luminosity in the form of non-thermal particles is 
-O.lxLj. 

5.2.2. Radiation 

Leptonic emission 

We can compute the SED of the synchrotron and IC com- 
ponents adopting a homogeneous model for the emitting re- 
gion. Since the magnetic field cannot be determined, two B- 
values, rendering different but representative situations, have 
been adopted: B = 0.4 and 50 G. The first value yields a dom- 
inant very high energy IC component reaching energies well 
above 1 TeV. The second value corresponds to a high magnetic 
field being in the limits of the hydrodynamical approximation. 
This latter case implies dominant synchrotron radiation. To com- 
pute the impact of synchrotron self-absorption, Dj s has been 
adopted as the width of the emitting region. 

For B = 50 G, U, ~ 300 erg cirT 3 , and 77 ~ 2 x 10~ 2 , 
synchrotron energy losses limit ~ 1 TeV, and the energy 
in non-thermal particles are radiated with similar luminosities 
via synchrotron and IC channels due to comparable timescales. 
Above the injection energy, the (cooled) particle spectrum fol- 
lows a power-law of exponent -2 up to the energy at which con- 
vection and radiative timescales become equal, around 1 GeV. 
In the one-zone approximation, electrons below 1 GeV escape 
the emitter before cooling completely, although they could still 
produce radiation outside the binary system. Above this energy, 
the particle spectrum becomes radiatively cooled with index -3. 



We have also explored the case of B — 0.4 G, for which the 
maximum energy reaches around 10 TeV. In this case, low mag- 
netic fields render a synchrotron component with luminosities 
well below the IC ones and a hard particle spectrum at the high- 
est energies due to the impact of KN IC losses, which harden 
the particle energy distribution when dominant (on the contrary, 
synchrotron and Thomson IC losses steepen the spectrum; see 
Khangulyan et al. 2008 for a thorough discussion on this). We 
show in Fig.[5]the corresponding SED for these two cases. The 
total radiated luminosity is ~ 10 35 erg s _1 , being shared by the 
synchrotron and the IC components with their proportion de- 
pending on B. It is worth noting that the age of the source pro- 
ducing the SED presented in Fig.|5]has been fixed to the typical 
lifetime of the shocks (~ R rb/Vj). 

We note that for B = 50 G, when synchrotron emission is 
high, the synchrotron self-absorption frequency is ~ 60 GHz 
when assuming a low energy cutoff at ~ m e c 2 , effectively yield- 
ing little radiation in the radio band. If the cutoff were at higher 
energies, electrons would radiate with a monoenergetic particle 
distribution in the radio band for a such magnetic field, rendering 
a very hard SED (of exponent 4/3) and little radiation in the ra- 
dio band again. Therefore, it is unclear whether radio emission in 
this scenario is detectable coming from the accelerator/emitter. 
Electron cooling and transport with the flow bulk may yield ex- 
tended synchrotron radiation in the radio band in case radio emit- 
ting electrons were carried to farther regions of appropriate B. 

IC scattering and photon-photon absorption in the stellar ra- 
diation field have been considered using angular averaged cross 
sections. Since we do not apply this model to any specific object 
at present, we do not account for the geometry of these interac- 
tions. Accounting for the geometry of photon-photon absorption, 
the latter can have less impact for certain accelerator/emitter lo- 
cated locations (k, 10 12 cm) and orbital phases. This would im- 
prove the detection chance for this kind of sources in the TeV 
range. We refer to Khangulyan et al. (2008) for an exhaustive 
discussion of IC and photon-photon absorption angular effects. 
Secondary pairs produced in the system may also be relevant 
from the radiative point of view (see Khangulyan et al. 2008; 
Bosch-Ramon, Khangulyan & Aharonian 2008b. 

Hadronic emission 

In case of protons, the densities in the jet are in general low 
and proton-proton energy losse^l are not very important, yield- 
ing little amount of 7r°-decay gamma-rays. The most suitable 
case is when proton escape velocity is low and density high in 
the region. This could occur in the shell region. Recollimation 
shocks would lead to jet density enhancement and low convec- 
tion velocities. Also, the mixing of jet and stellar wind material 
may also lead to high density regions that could suffer the impact 
of relativistic protons. 

In our context, gamma-ray fluxes <: 10 32 erg s would be 
expected from proton-proton interactions if « t /Vd ~ 10 (i.e. 
V A ~ 10 9 cm s~' implies n, > 10 10 cirT 3 ) for L p = 10 35 erg s _1 , 
where L p and n t are the relativistic proton luminosity and the 
target density, respectively. Although the efficiency of hadronic 
processes is relatively low for the production of gamma-rays, 
the interaction of powerful jets and dense winds could indeed 
lead to the generation of detectable TeV emission. In these 
proton-proton interactions, neutrinos (e.g. Romero et al. 2005; 
Aharonian et al. 2006a; Christiansen, Orellana & Romero 2006) 
and secondary electron-positron pairs would also be produced, 



5 Other hadronic channels, like synchrotron proton, photo- 
disintegration, or photo-meson production, have little efficiency 
or an extremely high threshold energy. They are not considered here. 
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Fig. 5. Typical computed SED of the synchrotron and IC emis- 
sion produced in a shock region. We adopt B = 0.4 (solid line) 
and 50 G (long-dashed line). The production (thin line) and the 
photon-photon absorbed (thick line) SED cases are presented. 



with luminosities similar to the gamma-ray ones. The radiation 
of secondary pairs would be relevant mainly at low energies (e.g. 
Orellana et al. |2007j l. 
Variability 

Initially, the forward and reverse shocks can produce non- 
thermal radiating particles as long as the jet goes through the 
system, with timescales ~ /? rb/^j ~ 300 s. In case recollima- 
tion shocks in the jet due to cocoon material pressure formed, 
particle acceleration would last as long as the cocoon pressure 
were higher than that of the jet. Since the cocoon pressure goes 
roughly like oc 1 ft 1 once the jet head leaves the binary system, 
cocoon recollimation shocks inside the system should not last 
longer than ~ R or b/Vj. All these events would appear as short 
bursts of X-ray and even GeV-TeV emission, and perhaps also 
in the radio band. 

Irregular jet ejection on timescales > R OI b/V w (the wind re- 
plenishment timescale) in a massive system could lead to (re- 
current?) very high energy flares as that observed in Cygnus X- 
1 (Albert et al. 2007) or apparently in certain orbital phases in 
LS 5039 (Aharonian et al. |2006bt and LS I +61 303 (Rico |2007l >. 
The observed ~ hour scale variability associated to these flares 
could be explained by a new ejection facing stellar wind replen- 
ished during a previous jet inactivity period of similar or longer 
duration. Eventual recollimation shocks due to wind ram pres- 
sure instead of cocoon pressure could stand as long as jet in- 
jection lasted and wind properties were relatively smooth. The 
variability of emission in this phase would be likely dominated 
by jet injection fluctuations or wind inhomogeneities, with asso- 
ciated timescales <: Rj/c. 

5.3. The fate of the jet 

As a test for jet bending, numerical simulations in slab geom- 
etry have been performed. It has been shown that the wind 
can efficiently transport shell and cocoon material sidewards. 
In the same way, the wind can significantly deflect the super- 
sonic jet itself via transferring of enough lateral momentum. We 
have shown that even in the case of a jet with kinetic power 



10 erg s , the wind seems able to disrupt the flow. Such a dis- 
rupted jet may look like a radio structure changing its morphol- 
ogy along the orbit, like that found by Dhawan et al. 12006 in 
LS I +61 303 and interpreted by these authors as the collision 
between the star and a pulsar winds. This interpretation would 
be challenged by the morphologies obtained from hydrodynam- 
ical simulations by Romero et al. (2007b) and Bogovalov et al. 
(120071) . 

The simulations and calculations presented here show that 
the recollimation shocks can, on one hand, easily develop inside 
the binary system region and, on the other hand, may influence 
significantly the evolution of the jet. As already stated, these 
shocks may generate efficient radiation and also could favor jet 
mass loading and deceleration via pinching, thus making deflec- 
tion by the wind easier. These shocks were already reported in a 
work by Peter & Eichler dl9951 l, in which the collimation of jets 
by the inertia of the ambient medium was studied. In the maps 
shown in that work, shocks similar to that observed in Fig.|4]are 
observed. 

The jet parameters adopted here result in very small growth 
rates of Kelvin-Helmholtz instability, mainly in the cases II and 
III. Therefore, disruption due to the growth of such an instability 
inside the region of interaction with the wind is not expected. 
Otherwise, our simulations show that the jets develop non-linear 
structures due to their initial overpressure, thus the growth of 
linear instabilities is not expected either. In short, the stability of 
the jet depends mainly on its inertia as compared to that of the 
wind and on the strength of the recollimation shocks (Perucho et 
al. 120051 [20071 . 



6. Summary 

The crossing of a hydrodynamical supersonic jet through the 
inner parts of an HMMQ renders strong "forward/reverse" and 
recollimation shocks. These shocks produce strong pressure and 
density enhancements, and some jet deceleration. The impact of 
the wind ram pressure on weak jets can prevent the formation 
of larger scale jets because of strong bending and jet disruption. 
We estimate that this occurs for jet kinetic luminosities below 
10 36 erg s _1 for the strong winds of O type stars. 

If efficient leptonic and/or hadronic particle acceleration 
takes place, the presence of photon, matter and magnetic fields, 
will yield significant amounts of radiation from radio to very 
high energies. On one hand, leptons would radiate mainly via 
synchrotron and IC processes. On the other hand, protons could 
interact with shell nuclei producing gamma-rays via neutral pion 
decay, and electron-positron pairs and neutrinos via decay of 
charged pions. Secondary pairs may also yield some amount of 
low energy radiation. In addition, the complexity of the wind/jet 
interactions leads to a plethora of variability timescales of the 
emission that could be more relevant than any intrinsic jet fluc- 
tuation. 

Further three-dimensional numerical simulations including 
more realistic implementations of the problem, magnetic fields, 
and radiative processes, should improve our understanding of the 
physics of jets in microquasars. 
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